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Abstract 

An integral formulation for acoustic radiation in moving flows is pre¬ 
sented. It is based on a potential formulation for acoustic radiation on weakly 
non-uniform subsonic mean flows. This work is motivated by the absence of 
suitable kernels for wave propagation on non-uniform flow. The integral solu¬ 
tion is formulated using a Green’s function obtained by combining the Taylor 
and Lorentz transformations. Although most conventional approaches based 
on either transform solve the Helmholtz problem in a transformed domain, 
the current Green’s function and associated integral equation are derived in 
the physical space. A dimensional error analysis is developed to identify the 
limitations of the current formulation. Numerical applications are performed 
to assess the accuracy of the integral solution. It is tested as a means of ex¬ 
trapolating a numerical solution available on the outer boundary of a domain 
to the far held, and as a means of solving scattering problems by rigid surfaces 
in non-uniform hows. The results show that the error associated with the 
physical model deteriorates with increasing frequency and mean how Mach 
number. However, the error is generated only in the domain where mean how 
non-uniformities are signihcant and is constant in regions where the how is 
uniform. 
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1. Introduction 


Predicting noise radiation from complex sources in moving flows is rele¬ 
vant to the automotive, energy and aeronautical industries. Noise radiation 
from turbofan nacelles and from other aircraft sources is a problem of particu¬ 
lar interest in the aviation sector. Numerical simulation of noise radiation and 
scattering can signihcantly reduce costs for design and certihcation. However, 
an efficient numerical method for high frequency noise propagation on non- 
uniform moving flows has not yet been demonstrated. Solving high frequency 
short wavelength problems on moving flows remains computationally expen¬ 
sive. In the aeronautical industry, noise propagation on non-uniform flows 
is typically predicted using finite element methods (FEM) P, discontinuous 
Galerkin methods (DGM) [2] and high order finite difference schemes [3]. 

Although volume based methods, such as FEM, DGM and finite difference 
schemes are able to solve wave propagation on a non-uniform flow, predict¬ 
ing noise radiation in unbounded domain requires the computational domain 
to be truncated. The truncation of the domain allows acoustic waves to be 
damped in a non-physical absorbing zone [31 111 E] and satisfy the radiation 
condition at the outer boundary of the domain. Moreover, these methods 
suffer of dispersion error and pollution effects [B]. These features are rel¬ 
evant limitations in case of noise radiation for large-scale short-wavelength 
problems. 

On the other hand, numerical methods based on boundary integral for¬ 
mulations, such as the boundary element method (BEM) [7] , inherently sat¬ 
isfy the radiation condition in the kernel and allow wave propagation in 
unbounded domains to be solved more effectively than in the case of vol¬ 
ume based methods. Moreover, the fast multiple BEM (FMBEM) is an effi¬ 
cient algorithm to solve wave radiation and scattering for large-scale short- 
wavelength problems [8]. However, BEM can only solve wave propagation 
exactly on uniform mean flows. Extending this method to non-uniform flow 
regions would be benehcial to a number of applications, such as forward fan 
noise acoustic installation effects. A surface integral formulation including 
non-uniform flow effects would also extend the applicability of wave extrap¬ 
olation methods. These approaches use an integral formulation defined on a 
closed surface on which the acoustic field is sampled from an ‘inner’ domain 
to radiate the solution to the far field. At the moment these methods are 
limited to uniform flow PHQ]. 

Gurrent boundary element modelling practices use the Lorentz transfor- 
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mat ion mm to solve wave propagation on a uniform mean flow. This 
variable transformation allows the uniform flow Helmholtz equation to be 
reduced to the standard Helmholtz problem without approximations. By 
means of a Lorentz transformation, BEM solvers for the standard Helmholtz 
equation can therefore be used for wave propagation on uniform flows. How¬ 
ever, due to the variable transformation, the physical space is deformed in 
the direction of the mean flow. The deformation of the domain complicates 
the formulation of the boundary conditions and the implementation of the 
transmission conditions for coupled formulations (nidi]. Alternatively, this 
drawback can be overcome by using an integral formulation in the physical 
space as proposed by Wu and Lee [U] in the frequency domain and by Hu 
PI in the time domain. 

For BEM, only approximate formulations are available for representing 
non-uniform mean flow effects. Astley and Bain im provided an approximate 
formulation for wave propagation on low Mach number mean flows based on 
Taylor’s transformation [TSl ITU] . In the same work, Astley and Bain PI re¬ 
ported an error analysis for Taylor’s wave equation showing that the accuracy 
of the physical model depends only upon the mean flow Mach number and 
the characteristic length scales of the acoustic waves and the mean flow. On 
the other hand, Tinetti and Dunn [20] provided a generalized local Lorentz 
transformation to represent the effect of non-uniform mean flows on wave 
propagation. However, the method has been restricted to mean flow helds 
with small gradients. Another approach is to move the terms including non- 
uniform flow effects to the right hand side of the equation and treat them 
as sources in the domain. The dual-reciprocity method (DRM) |2l] is then 
used to convert the domain integrals into boundary integrals. The absence 
of a robust method to dehne interpolating source functions for the DRM 
restricts the applicability of this approach. Thereby, modeling non-uniform 
flow effects for BEM is still an open problem. 

In this article, we present, in the physical space, an integral formulation 
with non-uniform flow based on a combination of the physical models asso¬ 
ciated with the Taylor and Lorentz transformations. The proposed physical 
model is an approximate formulation of the full linearized potential wave 
equation for isentropic compressible flows. The integral formulation derived 
applies to sound radiation on a weakly non-uniform potential mean flow. 
Consider the mean flow as a sum of a uniform and a non-uniform compo¬ 
nent which vanishes at inhnity. The term weakly non-uniform indicates that 
the non-uniform portion of the mean flow is small compared to the uniform 
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Q S ^ streamlines 



Figure 1: Schematic diagram of the reference problem showing mean flow streamlines in 
the solution domain 

part. A free field Green’s function is also determined for a subsonic weakly 
non-uniform flow as a kernel for the integral equation. Moreover, an error 
analysis is presented to extend and to revisit the estimate provided by Astley 
and Bain mi. This analysis shows the dependency of the error related to the 
physical model on the mean flow Mach number and on the frequency. The 
proposed formulation will be shown to improve the accuracy of the model 
compared to existing formulations based on either Taylor or Lorentz trans¬ 
forms being applied separately. 

The paper is structured as follows. Section presents the physical mod¬ 
els. In Section integral formulations are derived for wave propagation on 
weakly non-uniform mean flows and the Taylor formulation in the physical 
space. The Green’s functions associated with the integral formulations are 
derived in Section Boundary element formulations consistent with the 
proposed integral solutions are presented in Section]^ for an arbitrary source 
distribution. In Section a dimensional error analysis is developed to de¬ 
scribe the limitation of the proposed solutions. Finally, in Section some 
numerical results are presented to benchmark the integral formulations. 

2. Physical model 

A numerical solution to external noise radiation and scattering of a sound 
held from a source S' in a domain G by a body dQ in a non-uniform potential 
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subsonic mean flow is sought (see Fig. [^. Consider an inviscid, adiabatic 
and irrotational flow and assume that acoustic perturbations are of small 
amplitude compared with the steady mean flow. Under these hypotheses, 
wave propagation on a non-uniform mean flow can be described by means of 
a potential formulation |22] as, 


-Pq / Po 
Dt Dt j 


- V ■ 


0 


( 1 ) 


where 0 is the acoustic velocity potential and Do/Dt = d/dt + UQ-V denotes 
the material derivative over the mean flow; po is the mean flow density, cq the 
speed of sound and uq the mean flow velocity. Consider Bernoulli’s equation 

cS = cL-^(l|u„f-||u„f) (2) 

where Coo and Uoo are the speed of sound and the mean flow velocity vector 
in the far held. By means of Eq. ([^ and introducing the state equation, 

dpo _ 7 - 1 del 

n 9 

Po ^ Cq 

Eq. ([^ reduces to: 

^+2uo ■ + uo ■ V(uo ■ V0) 

+ ^V0 • V(uo ■ Uo) + (7 - 1) ^V ■ Uo (4) 

+ - i|uooinV^(^ = 0. 

Equation ([^ is rewritten as Eq. (|^ to allow a dimensional analysis to be 
performed. Mayoral and Papamoschou [23], following Astley and Bain im. 
have provided a description of the dependency of the terms in Eq. (|^ on the 
mean how Mach number M^o = Uoo/coo, the acoustic characteristic length 
scale La and the characteristic length scale of the mean how Lm- On this 
basis, Eq. (|^ is simplihed retaining only terms of order [(j)]/L\, Moo[(p]/L\, 
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and M^[0 ]/LaLm to give 


5^0 

dt'^ 


— d(j) 

.2u„.V^ 




-V0 ■ V(uo ■ Uo) + 


+ Uo ■ V(uo ■ V0) 

^ - 1 


|Uo| 


u„ 


2)V20 = O. 


(5) 


Assume that La < Lm and M^o -C 1. Then consider the mean flow formed 
by a uniform component Ugo and a non-uniform portion Ug. The main idea 
is to retain only first order effects in Ug over wave convection due to Ugo- If 
Uo = Uoo + Ug, where ||uq|| ||uoo|| and if it is assumed that the uniform 

mean flow velocity Uoo is aligned with the positive a:-axis, Eq. ([^ reduces to: 


<9^0 , ^ 2 V72I , 2 9^^ _n 

g^2 + 2Uo ■ V V 0 + Moo O' 


( 6 ) 


Hereafter, Eq. @ is referred to as the wave equation for weakly non-uniform 
potential flow. 

In the case of a uniform flow, where Uo = Uqo at all points in the flow, 
Eq. (|^ reduces to the uniform flow wave equation without approximation, 

Ir+Moo^^ 0-cLv20 = O. (7) 


dt 


dx 


If the mean flow is uniform, solutions of Eq. Q are also solution of Eq. Q 
without approximation. On the other hand, if terms of the order are 
neglected in Eq. (§ , it reduces to 

+ 2u„ ■ = 0. (8) 

In all that follows. Equation ([^ is referred to as the Taylor wave equation^ 
because it is consistent with the solution based on the Taylor transforma¬ 
tion dZllIH]. It is an approximation of order M^o of Eq. Q accounting for a 
non-uniform mean flow. 

If a time harmonic problem in which all perturbed quantities varies as e^^^* 
is considered, the acoustic velocity potential in the frequency domain can be 
defined as 0 = where 0 is a complex amplitude. In the Fourier domain 
Eqs. (§, 0 and ([^ are respectively referred to as the weakly non-uniform 
potential flow Helmholtz equation, the uniform flow Helmholtz equation and 
the Taylor-Helmholtz equation. 


9^0 
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3. Boundary integral formulation 

3.1. Weakly non-uniform potential flow Helmholtz equation 

By assuming a harmonic time dependence for 0, Eq. (j^ can be rewritten 

ecf - 2ifcMo ■ V0 + W - = 0, (9) 


where Mq is the mean flow Mach number vector, k = ujfcoo, and M^o is 
the uniform flow Mach number, provided that Moo = (d<foo;0,0). First, the 
reverse flow operator associated with Eq. ([^ is deflned considering a mean 
flow in the opposite direction of the actual flow fleld to give 

A;20 + 2zfcMo-+ = 0. (10) 


Following the general approach of Wu and Lee [15], a Green’s function G 


for the fundamental reverse flow operator Eq. (10) satisfies 


PG + 2iA;Mo ■ VG + 




x), 


( 11 ) 


where Xp is an arbitrary fleld point, x denotes the point source location and 
6 is the Dirac delta function. The Green’s function G represents the effect 
of a point source in a non-uniform flow whose direction is opposite to the 
actual base flow. G therefore differs from the Green’s function of the direct 
fundamental operator and is derived in Section 

To obtain an integral formulation for the weakly non-uniform potential 
flow Helmholtz equation, Eq. @ is multiplied by G and Eq. 0 by 0. 
Subtracting these two equations and integrating over the domain G yields 


In 


k^G + 2iA:Mo • VG + V"G - M, 
Glk'^cp- 2ifcMo ■ V0 + V V - M- 






oo 


dx"^ 


dV 


( 12 ) 


(j)S{xp — x)dV. 


The divergence theorem is applied to the l.h.s. of Eq. (12) in order to obtain 
an integral on the boundary surface of the domain dQ. Gonsistent with the 
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assumptions already made in the derivation of Eq. (§. some terms of higher 
order in M^c will be neglected. Consider the linear term in Mq in Eq. (12), 


Since [IT] 


/ (2ifc0Mo • VG + 2ikGyiQ ■ VcfydV 
Jn 

= [ 2ifc[V ■ (MoG0) - G0V ■ Mo]dE. 
Jn 

|Mo ■ V(Mo ■ Mo) 


(13) 


V-Mo = 


(14) 


0 1_7^(M2-M^)’ 

where Mq = ||Mo||, is of the same order as and since high order terms 
in Moo have already been neglected in Eq. (§ , the second term on the r.h.s 
of Eq. (13) can be dropped, giving 

[ 2ik[V ■ (MoG0) - G0V ■ Mo]dV ~ [ 2iA:Mo ■ nGcj)dS, (15) 
Jn Jan 

where n is the normal unit vector to d^l (see Fig. [^. Hence, Eq. (12) can be 
rewritten as 


0(Xp) = 


'an L 


dn ^ dn 


dS 


'an L 


2i/cMo ■ nG(j) + ( G— - (j)— ] 


dG 


(16) 


dS, 


where Ux is the x—component of the normal vector to the boundary surface. 
As shown by Wu and Lee [T3] the contribution of the boundary integral 
at inhnity is zero provided that 0(x) and Go(xp,x) satisfy the Sommerfeld 
radiation condition with flow. 

Equation (16) applies to any point Xp internal to the domain H. On the 
other hand, for any point on the boundary surface dQ, a limit approach |7| 
should be used to overcome the singularity of the integral formulation, which 
is integrable in the sense of Cauchy’s principal value. The following equation 


extends the uniform flow solution of Wu and Lee[l5] and is derived in Ap- 
pendix ~K 


G(xp)(/)(xp) = 


'an 


Ian 


dn ^ dn 


dS 


2ifcMo ■ nG0 + M^ ( G^ — ) Ha 


dx 


dx J 


(17) 


dS, 

















where 


(7(x,) = 


r f dGo 

J dQ \ On 


Ml^n,)dS 


Xp G 

Xp G d^l 


(18) 


and Go( 
defined in [Appendix A 


’-P5 ■ 


is the Green’s fnnction associated with the static operator 


In snmmary, Eq. (17) provides a bonndary integral formnlation for wave 
propagation on weakly non-nniform potential mean flows. It models first 
order non-nniform mean flow effects on wave propagation. The accnracy of 
the formulation depends on the deviation of the non-uniform flow Mq on the 
uniform component Mgc- Equation ([T^ is exact for wave propagation over 
a uniform mean flow and equivalent to the formulation presented by Wu and 
Lee [15] for a uniform flow. It reduces to the standard Helmholtz integral 
equation as Mq —)■ 0. 


3.2. Taylor-Helmholtz equation 

For the Taylor-Helmholtz equation, 

k‘^(j) — 2i/cMo ■ V0 -|- = 0, (19) 

an integral solution in the frequency domain is sought. The Green’s function 
for the reverse flow operator associated with the above equation is denoted 
by Gt, where 


k'^Gr + 2ifcMo ■ VGt + V^Gt = -5(xp - x) 


( 20 ) 


Gt differs from G and is derived subsequently in Section using a different 
set of assumptions consistent with a solution in the Taylor transformed space. 
Equation (19) is obtained from Eq. (|^ assuming 1 and La < Lm- 


Hence, with the same assumptions, the integral solution to Eq. (19) is derived 


from Eq. (16) by dropping the terms of order giving 


0(Xp) = 


I on 


^ d(p 

Gttt, - 

on 


dG‘ 


— 2iA)Mo ■ nG(r0 I dS. 
on J 


( 21 ) 


Equation (21) is written for a generic point internal to the domain Lt. How¬ 


ever, the integral becomes singular on the boundary surface dVt. On dVt 
the singularity is integrable in the sense of the Gauchy’s principal value. 
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Following the procedure showed in Appendix A a limit approach to the 
boundary surface is performed. In this case, the static operator associated 
to Eq. (20) is the Laplacian, V^Gqt = — <5(xp — x). The static operator is 


the same as for the standard Helmholtz problem. Therefore, the Cauchy 
principal value integral, obtained in this case, is the same as for the standard 
Helmholtz equation [Tj. By dehning 


C(Xp) = 


I, 


dGo^ 
dQ dn 


dS 


Xp G H 
Xp e d^l 


( 22 ) 


where Gqt = l/(47ri?) and R is the distance between the source and the 
observer, the following integral solution is obtained: 


C'(Xp)0(Xp) = 


'dn 


Utt; - 

on 


dG 


^ - 2ifcMo ■ nGrA dS. (23) 

on ) 


4. Green’s function for weakly non-uniform mean flows 

The Green’s function of the reverse flow operator, Eq. 0. is derived by 
means of a variable transformation. To determine G, a Lorentz transforma¬ 
tion m is applied following a Taylor transformation [18] of the fundamental 
reverse flow problem. 


d^G 

dR 




y^G + u 


d^G 

dx'^ 




(24) 


where G = Ge’^h In the Taylor-Lorentz space, Eq. (24) reduces to the stan¬ 


dard Helmholtz problem. Since the fundamental solution of the Helmholtz 
operator is well-known [7|, the Green’s function in the physical space, G, 
is retrieved by applying the inverse Taylor-Lorentz transformation to the 
Green’s function in the transformed domain. 

First, a Taylor transformation is applied including the non-uniform flow 


component Uq to Eq. (24). The independent variables are transformed as 


X = X, Y = y, Z = z, T = t 


4>'n(x) 


(26) 


where 4’o(x) is the mean flow velocity potential corresponding to the non- 
uniform flow part Uq, X = {X,Y,Z) and T denote the Taylor space-time. 
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whereas (x, y, z, t) denote the physical space-time. The differential operators 
in the Taylor space are then given using the chain rule as 


where Mq = Uq/coo- 
the Taylor space: 


— = — v = Vx-—— 
dt dT' ^ Coo dT' 


(26) 


Using the above equation, Eq. (24) can be rewritten in 


d^G 2 

Q'J'2 




dG 

W 


- 2CooMoo 


A 

Wr 


VxG 


M'^dG 

Coo dT 


M'oil 2 02(7 


-f c2 M2 

' ^OO-^ -^CXD 


d^G M',, d^G 


0X2 


0X0T 


+ 


02(7 

C2 0T2 


1 0M',,0(7\ 

Coo dX dTj 


clS{X,T). 

(27) 


Note that the Taylor transform brings the effect of the linear terms in Mq = 
(Mg_^, Mq^^, Mq J) inside the standard wave operator and introduces addi¬ 
tional terms associated with the mean flow non-uniformities in those de¬ 
pending on Moo- In the above equation, consider that MooMq and 

Vx ■ Mg is of order Mq (see Eq. (14)). Hence, neglecting higher order terms 
in Mg consistent with Eq. 0 . Eq. (|27D can be approximated by 


02(7 


- 2CooMoo 


^ dG 2 V. 


02^ 

xG + clMl— = cl5{X,T). (28) 


Equation (28) is nothing but the convected wave equation in the Taylor space. 

Assuming a uniform mean flow aligned with the x-axis but opposite to the 
actual direction of the mean flow Mqo, apply a Lorentz transformation m 
The independent variables are transformed as. 




Y = Y, Z = Z, T = Tf3^- 


M^X 

Coo /^oo 


(29) 


where X = (X,U, Z), T denote the Taylor-Lorentz space-time and (3oo = 
\/l — M^. The differential operators can be written as: 


A 


1 0 

PoodX 


Qrp /^OO 

Moo d 

/dooCoo dT 


A 

A 

Wy 


A 

w 


^ ^ ^ (30) 

dZ 0Z’ 
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Equation (30) is applied to Eq. (28) giving 


d^G 


cL^iG = ciS{X,T). 


(31) 


The l.h.s. of the above equation is the standard wave operator in the Taylor- 
Lorentz space and is independent of the mean flow. 


Since a steady state problem is considered, Eq. (31) is rewritten in the 


frequency domain. The Lorentz transformation introduces a time contraction 
dehned by the factor /S^o which, in turn, becomes a frequency dilation. In 
particular, the angular frequency in the Taylor-Lorentz space is a; = oj/f^oo- 
On the other hand, the Taylor transformation introduces only a time delay. 


Therefore, Eq. (31) can be rewritten as. 


PG +V|G = -5(X), 


(32) 


where k = oj/cao- For free held boundary conditions, the solution of Eq. (32), 
given by the monopole source solution, is 


G(X) = 


1 e 


—\kR 




(33) 


where R = V+ Y^ + Z^. 


Equation (33) is reformulated in the physical space, where a harmonic 
solution G(x, t) = G(x)e^'^^ is sought. An equivalent harmonic solution in the 
Taylor-Lorentz space is given by G(X, T) = G(X)e'‘^^. Applying the inverse 
Taylor-Lorentz transformation to the independent variables in G(X, T) yields 


G(X)e‘"^' = G(x)exp 


1^ ( /5oot 

Poc. 




-jf3o 


A 


OO o 


where 


G(x) = G(x)e 


Mqqx I ^n(^) 


= G(x)e^‘"* (34) 


(35) 


on the basis of Eq. (33), i.e.. 


Equations (34) and (35) give an explicit expression for the Green’s function 


exp 


G(x) = 


—ik 


^/x2+l3^(y^+z2) MqqX ’I‘o(x) 


+ 


47r 


(36) 
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The above equation is extended to a generic source position 
writing 

g-ifccTM 

AtiRm 



{xs,ys,Zs), by 


(37) 


where um = [Rm + Moo (a; - Xs)]//3l^ + [<ho(x) - <ho(x5)]/coo is the gener¬ 
alized reverse flow phase radius, extending the deflnition given by Garrick 
and Watkins [21], and Rm = y/(a; - Xs)^ + [{y - Vsf + {z - is the 

amplitude radius. 

In the case of a uniform flow, Eq. (37) is equivalent to the solution pro¬ 
vided by Wu and Lee 


On the other hand, for the reverse flow Taylor 
Green’s function Gt, the terms depending on M^ are neglected in Eq. (37). 
It follows 


G 


TlX,X^j = 


iAhT/v 


AtiR 


(38) 


where crMr ^ -R+['I»o(x)-<I>i)(x5)]/Cco, R ^ + (y - + (2 - 

and $0 denotes the total mean flow velocity potential. 


5. Boundary element formulation 

Gonsider the weakly non-uniform potential flow Helmholtz equation, Eq. ([^, 
with a generic distribution of harmonic sources g{'x.) in the domain G 


k‘^cj) — 2i/cMo ■ V( 


+ Vy-M^U=9(x). 


(39) 


Following the same procedure adopted in Section 3.1, Eq. (11) is multiplied 


by (j) and Eq. (39) by G. The difference of these equations integrated over 


the domain fl is rewritten using the divergence theorem to give 



nG(j) 

dS, 


dS 


(40) 


where G = G(xp, x) and 0 = 0(x) unless stated otherwise, x G dQ and Xp 
is an arbitrary point either in G or on dQ. In the above equation Mq and n 
are given along the boundary surface dG (see Fig. [^. 
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Equation (40) can be solved numerically by means of a collocation BEM 

can be defined. The major down- 


Alternatively, a variational formulation 
side of a variational statement is the increase in computational cost due the 
double integration over the boundary surface. A collocation formulation is 
provided below. 


First, Eq. (40) is rewritten considering that 


d(j) d(j) dn d(j) dr d(p dr] d(j) d(j) d(p 

ox on Ox or Ox or] Ox On or or] 


(41) 


where r and r] denote the coordinates along the unit tangent vectors r and r] 
on dQ, such that the normal vector to the boundary is given hy n = rj x t. 


Substituting Eq. (41) into Eq. (40) yields: 


C(xp)(/)(xp) = / G(xp,x)5((x) dV 
Jo, 

dd) BC 

- 2iA;Mo ■ 

Bn Bn 


'do 


dS 


(42) 


Mi 


'dn 


( d(j) 

-—1 


B(l> B(j) \ BG 

G 

\ Bn Bt Br] I Bx 


n^dS. 


Equation (42) is convenient from a computational point of view because 


the tangential derivative can be expressed as a sum of the shape functions 
multiplied by the nodal values of (j) as shown by Wu and Lee [T3] . 

Secondly, the above equation is discretized introducing a polynomial ex¬ 
pansion of the acoustic velocity potential and its derivatives [T5] 


0(x) 

NDoF 

r=l 

50(x) 

Bn 

NDoF „ , 

r=l 

(43) 

B(j){x.) 

NDoF ^-\T { \ 

Brjrix) 

NDoF / \ 

1 BNy.(x) ^ 

(AA\ 

Br] 

dr] 

r=l 

Br 

2^ dr 

r=l 



where Ar(x) represents the r-th polynomial shape function and NDoF is the 
total number of degrees of freedom. Hence, the discrete system of equations 


associated with Eq. (42) can be written as 


K'^(^ + K“^ =F, 

Bn 


(45) 
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where 


Kf^ = 5imC{^p,i) + / 2ik-Mo-rvGN^dS 
Jan 


'an 


dG , ,2 dG 


Nm + MiG 


dNm , dNm 

a '^X rs 

OT or] 


Ur > dS 


(46) 


T^U _ 

^Ira — 


'an 


^[1 - MiniN^dS, 


(47) 


Fi = / GgdV, 

Jn 


(48) 


G(xpx) is given in Eq. (37) and Xp ^ denotes the p-th collocation point. The 
vectors of the nodal degrees of freedom associated with the aconstic velocity 
potential and the aconstic particle velocity are denoted respectively as (p and 
drp/dn. 

Similarly, for the integral formnlation associated with the Taylor trans¬ 


formation, Eq. (40) can be approximated by neglecting terms of order Mi. 


Hence, using the Green’s function Gt yields: 


G(xp)0(xp) = / Gt(xp,x)5((x) dV 
Jn 


dp 


'an \ dn 


- 2ifcMo ■ ivGtP 1 dS. 
on 


The above equation can then be discretized using Eq. (43), but the derivation 
is omitted for the sake of brevity. 


6. Error estimate 

The weakly non-uniform potential flow wave equation, Eq. ( |^ , is an ap¬ 
proximation of the full potential linearized wave equation Eq. Q accurate 
to the hrst order in Mq. It is an exact formulation only for wave propagation 
on a uniform mean flow. This section presents a dimensional error analysis 
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of the weakly non-uniform potential flow wave equation compared to the full 
potential linearized wave equation. If La and Lm are respectively the char¬ 
acteristic length scales associated with the acoustic held and with the mean 
how held, Eq. Q can be expressed in the frequency domain as follows: 


u: 


— 2 ia;uo ■ Vd) -|- lLV 6 — u 


2 v 72 i „,2 




= E{ct>), 


(50) 


where 
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(51) 


In Eq. ([^ the terms on the l.h.s. of Eq. (50) are retained, whereas the terms 
included in the error function E are dropped by assuming Mq M^o- 

Following Astley and Bain na, Eq. is divided by and the r.h.s. 
terms which form E{4>) are rewritten as indicated below: 
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(62) 


where [ 0 ] is the characteristic dimension of 0 . 

To correctly bound the error brought by the solution of Eq. (17), a com¬ 
bined Taylor-Lorentz transformation is applied to Eq. (50). In the trans¬ 
formed space, 0 and G are exact solutions of the Helmholtz problem and the 
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fundamental Helmholtz operator. Following the procedure described in Sec¬ 
tion]^ to obtain Eq. (31), but considering the actual flow direction, Eq. (50) 
is transformed as follows: 


u 




= E, 


(53) 


where 


E =E -vlq- Uq— 0 + icu(Vx ■ Uo)0 - 2—Uoo ■ u^c; 


lOJ- 




dX 


+ OJ 


2 ^OjX^oo 




(54) 


Note that the above equation is a hybrid formulation. The term E is given 
in the Taylor-Lorentz space, whereas E is provided in the physical space and 


all of the other terms on the r.h.s. of Eq. (54) are written in the Taylor 


space. However, Eq. (54) is still suitable to estimate the order of magnitude 


of each term individually. The error terms in Eq. (54) introduced by the 

.2 

-"CO ’ 


Taylor-Lorentz transform can also be written, dividing again for c^, as 
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^ dX LaLm 

Hence, the error associated with the weakly non-uniform potential flow wave 
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equation Eq. (54) can be written as 
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^ A<M=oW , ^ ^ 
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* r 2 

^A 


ixi« 

(56) 


where Ci, (72, ■., (7g are constants of order 1. 

From Eq. ( |56[ ), the error E scales with 1/L^ = //coo where / is the 
frequency. In terms of mean flow length scale, the error varies with 1/Lm- 
The error therefore decreases as the mean flow becomes more uniform. The 
error vanishes in a uniform flow (Mq = 0). In other words, the accuracy 
of the formulation deteriorates only for wave propagation in a non-uniform 
region. 


The error e on the solution 0 to Eq. (53) is obtained by convolving E 
with the Green’s function (7, i.e.. 


e(xp) = / (7(xp,x)E(x) dV{yi). 


(57) 


In the Taylor-Lorentz space G varies as where Rm is the amplitude 


radius in Eq. (37). Given the geometrical length scale D associated with the 


domain hi, the above equation can be rewritten including Eqs. (52) and (55). 


The error contribution provided by the second and third terms in Eq. (52) 


to Eq. (57) is given by 


£1 ~ [0]MoMo, 


D 

Lm LaRm 


(58) 


The error Si increases linearly with frequency. However, if LaRm is constant, 
i.e. if the amplitude radius is inversely proportional to the wavelength, the 
error becomes independent of frequency. 


This applies also to all the remaining terms in Eqs. (52) and (55) except 
for the error terms of order 1/L\. In this case, the last term of Eq. (52) 
generates an error, 62 , in Eq. (1^ that is given by 


£2 ~ [0]MoMooT) 




L\Rm 


(59) 
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Note that the above equation scales quadratically with frequency. The con¬ 
tribution of this term is signihcant only at high frequency. 


Equations (58) and (59) can also be rewritten introducing the geometrical 


and acoustic relative distances of the observer, i.e. Rm/D and Rm/La, as 


El = [0]MoMoo — 

J^M 


D f D 


R 


M 


R 


M 


£2 = [0]MoMo 


D 


R 


M 


R 


M 


(60) 

Small values of Rm /D denote the geometrical near field whereas large values 
of Rm/D are associated with the geometrical far held. Similar dehnitions can 
be given for the acoustic held on the basis of Rm/La- The above equations 
show that, when sound propagates in a non-uniform how, the error decreases 
in the geometrical far held but increases in the acoustical far held. Contours 
of El and £2 on a logarithmic scale are shown in Figs. and 
error analysis can be extended to 2D problems, where G scales with 1/^/kRM, 
but it is not explicitly reported here for the sake of conciseness. 


2b The present 




Rm/La [n.d.l Rm/La [n.d.l 


(a) £1 (dB) (b) £2 (dB) 

Figure 2: Contour of the error estimates, £i and £2, on the weakly non-uniform potential 
flow equation Eq. ([^ against the full potential linearized wave equation Eq. Q. 

The same conclusion holds for the Taylor-Helmholtz formulation. How¬ 
ever, the error increases even when wave propagation occurs on a uniform 
how. This is because the Taylor-Helmholtz equation is also an approximation 
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of the uniform flow convected Helmholtz equation, namely the error E does 
not vanish in a uniform flow region. 

7. Benchmark problem 

In this section, numerical examples are provided to assess the accuracy 
of the boundary integral solution to the weakly non-uniform potential flow 
Helmholtz equation Eq. ([^ compared to a reference solution of the full po¬ 
tential linearized Helmholtz equation based on Eq. Q. A comparison with 
integral solutions to the uniform flow Helmholtz and Taylor-Helmholtz equa¬ 
tions is also provided. 

First, the integral formulations are used to solve a ‘wave extrapolation’ 
problem where an ‘inner’ solution is obtained on an arbitrary closed surface 
in the flow and extrapolated to the far held by using the proposed integral 
formulation. Second, a traditional boundary element (BE) solution is pre¬ 
sented. The scattering of a sound held by a rigid body is computed at all 
points external to the scatterer up to the boundary surface. In this study, 
the numerical issue associated with irregular frequencies in the BE solution 
is overcome by avoiding the characteristic frequencies of the scatterer. 

The scattering by a rigid cylinder of radius a of the acoustic held generated 
from a monopole point source in a non-uniform mean how is used as the 
basis for both the test cases. A full 2D problem is solved. The analytical 
solution of a potential inviscid incompressible mean how around a cylinder 
without circulation is used to dehne the base how and a monopole point 
source of unit magnitude is dehned at a point S. The problem is solved in an 
unbounded domain for far held Mach numbers in the range 0.0 —)■ 0.3. The 
mean how density poo = 1-22 kg m“^ and the speed of sound c^o = 340 m 
s“^ are constant. The reference pressure for the computation of the SPL = 
20 log]^Q(Pfms/Pre/) is Pref 2 ■ 10 Pa. 

As a reference solution, a Lagrangian FEM with cubic element interpo¬ 
lation based on 30 degrees of freedom per wavelength is used to solve the 
full potential linearized wave equation Eq. ([^. The problem is solved in 
the frequency domain and a perfectly matched layer (PML) [5] is applied 
to model the radiation condition. The PML is located at r = 10a. Due to 
the approximation of the source model and the radiation condition, an error, 
measured against an analytical solution |26], of 0.05% at the outer surface of 
the computational domain is assessed to be provided by the EE solution in 
the case of quiescent media. 
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7.1. Wave extrapolation test case 

The first test case consists of a monopole point sonrce located at Xg = 
(0, —1.5a) as shown in Fig. In this case, the integral formnlation is nsed to 
extrapolate the solntion to the far held. First, the reference FE solntion of the 
fnll potential linearized Helmholtz eqnation solves radiation and scattering 
of a monopole sonrce by the cylinder with a non-nniform how in the inner 
domain Second, the reference solntion is sampled on a closed control 
snrface dflin, inclnding the cylinder and the sonrce, and radiated to the far 
held Qout by means of an integral formnlation. In the cnrrent instance dQin 
is located at a radial distance Vcs = 2a. The maximnm valne of Mq on 
dQin is Moo/4. The order of interpolation of the integral solntion on flin is 
consistent with the hnite element model. The FE reference solntion is also 
nsed as benchmark in the far held. 



Fignre 3: Geometry of reference for the wave extrapolation test case. Scattering of 
the sound field from a monopole source on a potential mean flow by a rigid cylinder. The 
sound field in the outer domain flout is extrapolated based on the field on the inner surface 
BQ. 

(_/a LiYi • 


First, the error generated by the integral solntion associated with the 
weakly non-nniform potential how Helmholtz eqnation against the reference 
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Figure 4: error for cj) along a field point circular arc with radius rfp = 8a a.s a function 

of the non-dimensional frequency ka for the wave extrapolation test case of Fig. The 
solution is based on the integral formulation for the weakly non-uniform potential flow 
Helmholtz equation Eq. (17). 



Figure 5: error for (f) along a field point circular arc with radius rfp = 8a as a function 

of the mean flow Mach number Moo for the wave extrapolation test case of Fig. The 
solution is based on the integral formulation for the weakly non-uniform potential flow 
Helmholtz equation Eq. (17). 
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Figure 6: error for (j) at ka = 4.25 along field point circular arcs with radii rfp = 6a, 8a 

and 10a as a function of the mean flow Mach number for the wave extrapolation test 
case of Fig. The solution is based on the integral formulation for the weakly non-uniform 
potential flow Helmholtz equation Eq. (17). 


FE solution is calculated. The error is defined as 


Ei2 = 100 X 


frjp ||0(x) -(^(x),e/Pd5 


\ IvfpUrefi^WdS 


(61) 


where Tfp denotes a closed circular arc of field points (see Fig. [^. The 
error for the acoustic velocity potential is shown in Fig. [^plotted against the 
Helmholtz number ka on a circular arc of field points with radius rfp = 8a. 
Note that the error is computed considering the physical distance of the 
observer R in hen of the amplitnde radius Rm since for Moo < 0.3 the 
maximum difference between R and Rm is abont 5% of R. The variation of 
the error for 0 against Mgo is shown in Fig. [^whereas the sensitivity of 

the error on the distance to the observer is shown in Fig. In the latter, 

the error is compnted on circular arcs of radii rfp equal to 6a, 8a and 10a for 
ka = 4.25. Note that the mean flow is almost nniform at these distances from 
the cylinder. To fnrther assess the accuracy of the integral formnlation for 
the weakly non-nniform potential flow Helmholtz equation, the compntation 
of the error for 0 against is also performed for the formnlation of Wu 


and Lee na and the Taylor-Helmholtz eqnation Eq. (H^. This comparison 
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is given in Figs. and [^computing the error along an arc of field points 

with radius rjp = 8a for different values of the the non-dimensional frequency 
ka. 





(a) ka = 1.94 


(b) ka = 4.25 


(c) ka = 9.24 


Figure 7: error for (j) along a field point circular arc with radius Vfp = 8a as a 

function of the mean flow Mach number for the wave extrapolation test case of Fig. 
at ka = 1.84, 4.25 and 9.24. The solutions are based on the weakly non-uniform potential 
flow Helmholtz integral formulation Eq. 0 (A), the uniform flow Helmholtz integral 
formulation [TS] (o) and the Taylor-Helmholtz integral formulation Eq. (231 (x). 


A description of the local accuracy provided by the three integral for¬ 
mulations is presented in Fig. which illustrates the pressure directivity at 
Tfp = 8a for ka = 9.24 and Moo = 0.3 plotted against the angle 6 ^ where 
6 ^ = 0 corresponds to the x-axis and 9 is measured counterclockwise. The 
acoustic pressure is given by 

p =-po(ic^0 + uo • V0). (62) 

Again, the FE solution of the full potential linearized Helmholtz equation is 
used as a reference result. The solution provided by the Taylor-Helmholtz 
equation introduces local errors of up to 10 dB. On the other hand, the 
error given by the weakly non-uniform potential flow Helmholtz solution and 
the related uniform flow approximation is limited to 5 dB. However, the 
weakly non-uniform solution improves the results given by the uniform flow 
Helmholtz equation in the shielded area [9 = 60° — 120°], where the mean 
flow is not aligned with the uniform stream. 

The accuracy of the proposed integral formulation Eq. 0 decreases 
linearly with frequency (see Fig. and Mach number M^o (see Fig. |^, 
whereas the error is almost constant when wave propagation occurs on a 
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Figure 8: Acoustic pressure directivity along a field point circular arc with radius rjp = 
8a for ka = 9.24 and M^o = 0.3. The wave extrapolation solution (see Fig. is based on 
the integral formulations Eq. ( [T7| ) (IF weakly non-uniform), Wu and Lee [T^ (IF uniform 
flow) and Eq. (231 (IF Taylor). 


uniform mean flow (see Fig. since the integral formulation associated with 
the weakly non-uniform potential flow Helmholtz equation is exact for a 
uniform base flow. These results validate the error analysis of Section 

Moreover, the integral formulation based on the weakly non-uniform po¬ 
tential flow Helmholtz equation provides the best accuracy independently 
of frequency and Mach number compared to the solution obtained using 
the integral formulations for the uniform flow Helmholtz and the Taylor- 
Helmhotlz equations (see Fig. [^. The non-uniform flow effects in the weakly 
non-uniform formulation improve the prediction made compared to the uni¬ 
form flow Helmholtz equation. The Taylor formulation however performs 
poorly and is less accurate than the uniform flow formulation which assumes 
uniform flows at all points. The problem with the Taylor equation is that 
the error grows even in the uniform flow region whereas the uniform and the 
weakly non-uniform flow models are exact once this domain is reached. 


1.2. Boundary element test case 

In this section, the accuracy of the boundary integral formulation for 
the weakly non-uniform potential flow Helmholtz equation is assessed by 
considering a direct BE solution of Eq. (42), solving for the surface acoustic 
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Figure 9: Reference domain for the BE solution showing the non-uniform mean flow 
streamlines. Scattering of the sound field from a monopole source S' by a rigid cylinder in 
an unbounded domain with non-uniform flow. 


velocity potential on the cylinder due to a sound field from a point source 
(see Fig. |^. The integral formulation is used to represent wave radiation 
and scattering up to the boundary surface, where the non-uniform mean 
flow component Mg is of the same order of magnitude as M^. The cylinder 
with radius a is centered at the origin of the reference frame and a monopole 
source is located at Xs = (—2a, —2a) as shown in Fig. The error in the 
solution is calculated against the reference FE solution. A cubic Lagrangian 
BE interpolation with 10 DoFs/A is used to solve the integral equations. 

The dependency of the error for 0 against the non-dimensional fre¬ 
quency at a circular arc of field points with radius r/p = 4a is illustrated 
in Fig. IT Again, note that the error is computed considering the physical 
distance of the observer, denoted previously as R, in lieu of Rm since for 
Moo ^ 0.3 the maximum difference between R and Rm is about 5% of R. 
Figure 11 shows the sensitivity of the error for </> to the distance of the 
observer sampling the error on circular arcs where the flow is almost uniform, 
for a non-dimensional frequency ka = 5. 

A comparison of the accuracy of the different integral solutions is also 
shown in Figs. 12a , 12b and 12c where the error for 0 is computed at 
an arc of field points with radius Vfp = Aa for the weakly non-uniform po¬ 
tential flow Helmholtz equation, the uniform flow Helmholtz equation na 
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Figure 10: error for (jj along a field point circular arc with radius rjp = 4a for the 

BE solution of the problem in Fig. as a function of the non-dimensional frequency ka. 
The solution is based on Eq. ([4^. 



Figure 11: error for at a non-dimensional frequency ka = 5 along field point circular 

arcs with radii rfp = 6a, 8a and 10a for the BE solution of the problem in Fig. as a 
function of the mean flow Mach number M^o- The solution is based on Eq. ([4^. 
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and the Taylor-Helmholtz equation. Figures 13a, 13b and 13c show contours 


of the real part of (j) over the solution domain for ka = 10 and M^o = 0.3. 


The weakly non-uniform potential flow formulation (Fig. 13a) clearly ap¬ 


proximates the full potential FE solution (Fig. 13b[ ) more accurately than 
the uniform flow Helmholtz solution (Fig. |13c ). This is more evident if the 
sound source is in a region where the mean flow is strongly non-uniform, as 
shown in Fig. where the same problem as in Fig. [^is solved locating the 
monopole point source at = (—1.3a, —0.5a). 



M,... [n.d.] 


|n.d.] 


(a) ka = 2 


(b) ka = 5 


(c) ka = 10 


Figure 12: error for (p along a field point circular arc with radius r/p = 4a as a 

function of the mean flow Mach number at non-dimensional frequencies fca = 2, 5 
and 10 for the problem in Fig. The solutions are based on the weakly non-uniform 
potential flow Helmholtz equation Eq. (42) (A), the uniform flow Helmholtz equation [T^ 
(o) and the Taylor-Helmholtz equation Eq. (49) (x). 


The previous results have been presented in terms of computed value of 
the acoustic velocity potential. The acoustic pressure p is another quantity 
of practical interest. The absolute value of the acoustic pressure at = 4a 
for Moo = 0.1 and 0.3 is shown in Figs. 15 and 16 The reference solution 


is compared to the approximate integral solutions at ka = 10 for different 
observer angular positions (the angle 6 measured from the x-axis and positive 
counterclockwise). For M^o = 0.1, all of the integral formulations are mostly 
within 1 dB and 5 dB of the reference solution. However, at M^o = 0.3 the 
error incurred in the Taylor-Helmholtz solution reaches 10 dB. Nonetheless, 
the BE prediction based on wave propagation on uniform flows overestimates 
the reference solution in the shielded area [6 = 30° — 70°], where the incident 
and the scattered held interfere destructively and the how is not aligned 
with the uniform stream. Note that the uniform how Helmholtz formulation 
neglects wave refraction ehects due to non-uniformities. 
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Figure 13: Real part of the acoustic velocity potential, Re{(j)), at a non-dimensional 
frequency ka = 10 for = 0.3. The BE solution of the weakly non-uniform potential 
flow Helmholtz equation (a), the FE solution of the full potential linearized Helmholtz 
equation (b) and the BE solution of the uniform flow Helmholtz equation (c) are 
shown for the problem in Fig. 



X [n.d.] 


r 5e-05 

- 4e-05 

- 3e-05 

- 2e-05 

- 6e-06 ^ 

t -6e-06 ^ 
-2e-05 
-3e-05 
-4e-05 
-5e-05 





5e-05 

4e-05 

3e-05 

2e-05 

6e-06 ^ 

-6e-06 

-2e-05 

-3e-05 

-4e-05 

-5e-05 



■p 5e-05 

- 4e-05 

- 3e-05 

- 2e-05 

- 6e-06 ^ 
-6e-06 
-2e-05 
-3e-05 
-4e-05 
-5e-05 




(a) 


(b) 


(c) 


Figure 14: Real part of the acoustic velocity potential, Re{(j)), at a non-dimensional 
frequency ka = 10 for M^a = 0.3. The BE solution of the weakly non-uniform potential 
flow Helmholtz equation (a), the FE solution of the full potential linearized Helmholtz 
equation (b) and the BE solution of the uniform flow Helmholtz equation [15] (c) are shown 
for the problem in Fig. |^but locating a monopole point source at = (—1.3a, —0.5a). 


For the weakly non-uniform potential flow Helmholtz equation the error 
varies linearly with frequency (see Fig. 10) and with (see Fig. 12). The 
error is independent of the observer distance when it is sampled on a region 
where the flow is almost uniform (see Fig. 0 because the formulation is 
exact for a uniform flow, namely no additional error is generated as the held 
point radius increases. These results validate further the error analysis of 
Section jO] 

As in the previous test case, the weakly non-uniform how Helmholtz 
formulation outperforms both the uniform how Helmholtz and the Taylor- 
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B [deg] 


Figure 15: Acoustic pressure directivity along a field point circular arc with radius 
r = 4a at a non-dimensional frequency ka = 10 for a mean flow Mach number = 0.1. 
The BE solutions (see Fig.|9) are based on the integral formulations Eq. ( [4^ (BE weakly 
non-uniform), Wu and Lee [IK] (BE uniform flow) and Eq. (49) (BE Taylor). 


Helmholtz solutions (see Figs. 12, 15 and 16). The improvement against 
the uniform flow Helmholtz solution is larger at low Mach numbers. Since 
M' Moo, for increasing M^o the error term of order M^Mq/L\ in Eq. ( [5^ , 
is signihcant both for the uniform and weakly non-uniform flow solutions, 
eventually reducing the improvement given by the non-uniform flow formu¬ 
lation on the uniform flow approximation. 


8. Concluding remarks 

A novel integral formulation for sound radiation on non-uniform mean 
flows has been proposed. It provides an approximate solution of the full 
potential linearized Helmholtz equation by assuming a weakly non-uniform 
potential mean flow. A Green’s function, describing non-uniform flow effects 
on wave propagation, was derived for potential subsonic mean flows. The key 
advantage is that the Green’s function and subsequent integral equation are 
presented in the physical space. The integral formulation is exact for wave 
propagation on a uniform flow and the error generated in the non-uniform 
flow region remains constant when wave propagation occurs on a uniform 
flow. This distinguishes the current approach from the Taylor transforma- 
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9 [deg] 


Figure 16: Acoustic pressure directivity along a field point circular arc with radius 
Tfp = 4a at a non-dimensional frequency ka = 10 for a mean flow Mach number = 0.3. 
The BE solutions (see Fig. [9 1 are based on the integral formulations Eq. (42) (BE weakly 
non-uniform), Wu and Lee [IK] (BE uniform flow) and Eq. (49) (BE Taylor). 


tion integral formulation applied by a number of previous researchers in 
computational Aeroacoustics (nmniEa]. 

A comparison with the Taylor integral formulation and the integral solu¬ 
tion for the uniform flow Helmholtz equation (Lorentz formulation) has been 
presented. First, an analysis based on wave extrapolation on a non-uniform 
flow was performed. In this case, the condition that Mq -C M^o was satis- 
hed at all points outside the integral surface. For the weakly non-uniform 
potential flow Helmholtz solution the error varies linearly with M^o and is 
proportional to the frequency /. The new approach performs better than 
the Taylor formulation in which the error varies as and is linear in /. It 
also outperforms the uniform flow Helmholtz equation, reducing the error of 
40% for Moo < 0.3. 

In a second test case the acoustic held was solved up to the solid scat¬ 
tering surface by means of a full BE approach. On the boundary surface 
itself Mq is of the same order of magnitude of Moo, pushing to the limit 
the underlying assumption of the formulation of the analysis Mq -C Moo- 
An error analysis on the integral solution based on the weakly non-uniform 
potential how Helmholtz equation shows that the error is proportional 


31 











to and varies linearly with /. A consistent improvement in accnracy 
is achieved against the integral solntions based on nniform mean flows and 
on the Taylor-Helmholtz eqnation for a non-nniform flow. A significant ad¬ 
vantage of the cnrrent weakly non-nniform potential flow Helmholtz solntion 
over the nniform flow Helmholtz solntion is evident at low Mach nnmbers, 
where for instance a redaction of 50% of the error is observed at = 0.1. 
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Appendix A. Integral formnlation along the boundary surface 


An extension of Eq. (|16|) to the bonndary snrface is given on the basis 

When an exterior wave 


of a limit approach [7], following Wn and Lee 
propagation problem is considered, the domain H (see Fig. is modified by 
snbtracting a hemisphere of radins e in the neighbonrhood of Xp G H, 


0(Xp) = 




/an+aoE . 


dS 

on on J 

2i/cMo ■ nG0 -|- I 


dx 


dx J 


(A.l) 


dS. 


When the radius e tends to zero, the surface of the hemisphere tends to zero 
as e^. Since G varies as 1/e, the contribution of the terms including G to the 
integral over the infinitesimal surface is zero. Hence, from Eq. (A.l) one has 

r [ . / f)rh nn\ i 

dS (A.2) 


'an 


- </(Xp) 


'aO( 


fdG ^.dG 


dx 

X I dS. 


dx J 
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In the limit of e —)■ 0 the Green’s function in Eq. (11) tends to the Green’s 
function of the static operator na Go- The static operator associated with 
Eq. 0 is given by 


V^Gn 


J,f2 « 


(A.3) 


The solution of the equation above for a free held problem gives Go(xp, x) = 
l/(47ri?M)- Wu and Lee [15] have shown that the integral on the hemisphere 
of radius e can be rewritten on the boundary surface dQ as follows: 


'dfi. 


OGq 

On 


j./r2 ^^0 


dS = 


'do. 


OGq 

dn 




dS. (A.4) 


Therefore, applying Eq. (|A.4|) to Eq. (|A.2|) yields 

G(xp)0(xp) = / 


I an 


'an L 


90 dG\ 

Gj. - dS 

on on J 

2iA;Mo ■ nG0 + 




dG\ 

j 


(A. 5) 


dS, 


where 


i 1 - hn (1“ - dS 


Xp G G 
Xp G 9G 


(A.6) 


Appendix B. 2D Green’s function 


This section presents the 2D Green’s function for wave propagation on 
a weakly non-uniform potential mean how. The formulation is limited to 
subsonic hows. The solution is based on the proof given in Section The 
Green’s function for the Taylor-Helmholtz formulation is also derived. 

The Green’s function for the Helmholtz operator in the Taylor-Lorentz 


space, Eq. (32), is given by 


G(X) = ' I kR.,o 


(B.l) 


where R 2 D = V -|- Y^, is the Hankel function of the second type of 
order zero and k = fc//3oo- A time harmonic solution is sought for G, namely 
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(^(X, T) = G'(X)e*‘^^. The inverse Taylor-Lorentz transformation is applied 
to the above equation retrieving the Green’s function in the physical space, 
as described in Section 


G(x) 


Vc " V j 


(B^2) 


Extending the Green’s function to a generic source position x^ = {xs,ys) 
yields 


G(x.x,) = —ff, 


( 2 ) 


kRM,2D\ 


(B.3) 


where Rm, 2 D = \/(x — XsY + ^'^{y — VsY- The Green’s function associated 
with the static operator Eq. (A.3) is Gq = —log{RM, 2 D)/‘^'^- 

For a uniform flow, Eq. (B.3) reduces to the expression provided by Bailly 
and Juve 123 including a reverse mean flow. On the other hand, if -C 1, 
the 2D Green’s function can be rewritten for the Taylor-Helmholtz formula¬ 
tion as 

G.(x.x,) = ^ffr(M,o)e-““^ (B.4) 

^Poo 

where R 2 D = y/(a; - Xs^ + (y - ysY- 
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